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The backreaction of inhomogeneities on the global expansion history of the Universe suggests a 
possible link of the formation of structures to the recent accelerated expansion. In this paper, the 
origin of this conjecture is illustrated and a model without Dark Energy that allows for a more 
explicit investigation of this link is discussed. Additionally to this conceptually interesting feature, 
the model leads to a ACDM-like distance-redshift relation that is consistent with SN data. 
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Averaged equations for the expansion of the 
Universe 

The averaging problem 

The evolution of our Universe may be described 
in terms of Einstein's equations of General Relativ- 
ity. These are ten coupled differential equations for 
the coefficients of the metric that describe our space- 
time. In the cosmological case, where we are only in- 
terested in the overall evolution and not in the detailed 
local form of the inhomogeneous metric, cosmologists 
widely work with the assumption that the global evo- 
lution is described by a member of the class of homoge- 
neous and isotropic solutions of Einstein's equations. A 
homogeneous-isotropic fluid is considered as source term 
in a Robertson Walker (RW) constant-curvature geome- 
try, determining the kinematical and dynamical proper- 
ties of the Universe in a scale-independent way. These 
properties are thereby condensed into the evolution of 
a single quantity, the scale factor a{t). The fundamen- 
tal question, dating back to [22j and most prominently 
raised by George Ellis in 1983 [lOj, is then, whether this 
procedure leads to the correct description of the global 
behavior of our spacetime. To address this question, 
one has to flnd a way to explicitly average Einstein's 
equations. This is necessary, because if one performs an 
average of an inhomogeneous metric whose time evolu- 
tion has been determined by the use of the ten Einstein 
equations, one flnds a result that in general differs from 
the classical model. This is obvious from the fact that 
Einstein's equations are nonlinear. But already at the 
linear level, as soon as the volume element of the do- 
main of averaging is time-dependent, this difference is 
present. This is easy to see from a derivation of the 
deflnition of the volume-average of a scalar quantity /, 
{f)v(t) ■■= Iv f{t,X)^(^)g{t,X)A^X/Vv(t) with re- 
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spect to time which provides 

dt if) J, - {dt f)r, = if - {0)j, 7^ , (1) 

being the local expansion rate. So, if we take an im- 
plicitly averaged metric, i.e. the RW metric, and then 
evolve a homogeneous matter source through the stan- 
dard Friedmann equation, this will not give the average 
of the time-evolved local quantity. Or, to put it short, 
time evolution and averaging do not commute [llj, of- 
ten strikingly written as G^y{{g^y)) / {G^y{g^j,)) and 
illustrated in Fig. 1. 



Provenance of the averaged equations 

In the recent literature, the question of how big the 
difference between these two approaches is, has received 
growing interest, mainly due to attempts to relate it 
to the volume acceleration (Dark Energy problem) [5j, 
[12],p!5( IT6],[6]. Since then, there have been many cal- 
culations trying to quantify the impact of this non- 
commutativity of time-evolution and averaging. The 
direct way of trying to average the Einstein equations 
in their tensor ial form is very involved, because it is 
not clear how to deflne a meaningful average of tensors. 
Therefore, a restriction to the averaging of scalar quan- 
tities is a way out of this problem and provides general 
features of the averaged dynamics [3','4]. In this approach, 
one uses the ADM equations to perform a 3 + 1 split of 
the spacetime into spatial hypersurfaces orthogonal to 
the fluid flow. The source is also in this case often taken 
to be a perfect fluid and the equations take their simplest 
form in the frame of an observer comoving with the fluid. 
After this split one can identify scalar, vector and tensor 
parts of the resulting equations. For the scalar sector 
there is then a straightforward deflnition of an average 
quantity as the Riemannian volume integral of the scalar 
function, taken over a mass-preserving comoving domain 
of the spatial hyper surface, and divided by the volume of 
this domain. The use of this deflnition astonishingly pro- 
vides a set of differential equations for the average scale 
factor of the averaging domain, that closely resembles the 
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FIG. 1: Illustration of the cause of the departure of the average evolution from the standard Friedmann solutions: the non- 
commutativity of spatial averaging and time-evolution. The standard Friedmannian picture is the left branch, where first the 
average over the inhomogeneous matter distribution is performed and the evolution of this homogeneous distribution is then 
evolved by the Einstein equations. In contrast to that, in the averaged model shown on the right branch, the evolution of the 
variables is calculated with the full metric and the average is taken over the evolved inhomogeneous distribution. The two 
approaches give in general different results. But the main question remains: How relevant is this difference? 



Friedmann equations in the homogeneous case: 
3Hl = SttG (g)^ - \ (n)^ -^Qv^A; 

{gj + 3Hv {g)j, = ; a^\alnv) + a^'(4Qi^) = 0. 

(2) 

This is surprising, because in this approach it is not nec- 
essary to constrain the matter source to a homogeneous- 
isotropic one, but one can have arbitrary spatial varia- 
tions in the density and the geometrical variables. There 
are, however, two important differences between the gen- 
eral averaged evolution equations for the volume scale 

1 /3 

factor ax> oc Vj^ and the Friedmann equations. First, 
there is one extra term Qx>, which is called the kinemat- 
ical backreaction term. It encodes the departure in the 
spatial hypersurface from a homogeneous distribution. 
It is defined as the variance of the local expansion rate 
of the spacetime minus the variance of the shear inside 
the domain V. Therefore, if the expansion fluctuations 
are bigger than the shear fluctuations, Qd is positive and 
contributes to the acceleration of the spatial domain. For 
dominating shear fluctuations, Qx) is negative and decel- 
erates the domain's expansion. This effective term Q^, 
that emerges from the explicit averaging procedure, also 
induces the second difference to the standard Friedmann 
equations: an integrability condition assures that the ex- 
pansion law is the integral of the acceleration law. While 
for the Friedmann equations this integrability is guaran- 
teed through mass conservation, here there is a gener- 
alized conservation law that, besides mass conservation. 



dynamically links Qx) to the average intrinsic curvature 
on the domain V. Unlike in the Friedmann case, where 
the curvature scales as CLx) 1 

the dependence of the aver- 
age curvature (7^)^ on ax) has not necessarily the form of 
a simple power law. In fact it can be shown that the cur- 
vature picks up a time-integrated contribution along the 
history of the evolution of Qx>^ and evolves in this way 
generically away from flat initial conditions, expected to 
emerge from inflation. 



Uncommon properties of the averaged model 

These two changes to the standard Friedmann equa- 
tions may alter the expansion history considerably. Re- 
garding Eqs. (2), it is easy to see that for Qx > 47rG {g)^^ 
there may even be an accelerated epoch of expansion 
without the presence of a cosmological constant. It may 
seem surprising that even in a Universe filled with a per- 
fect fluid of ordinary (or dark) matter only, there may 
be an effective acceleration of a spatial domain V. This 
occurs, because in the calculation of the average expan- 
sion rate of P, the deviations of the local expansion rates 
from the average introduce a positive-definite, non-local 
fluctuation term ((6> — {0)x)'^)x- Also, the average is 
weighted by its corresponding volume. Therefore, faster 
expanding subregions of P, which will by their faster 
growth occupy a larger and larger volume fraction of P, 
will eventually dominate the global expansion. Subre- 
gions whose expansion slows down will finally occupy a 
much smaller fraction of the volume of V. This implies 



3 



that a volume-weighted average of the expansion rate 
will start from a value between the one of the slow and 
fast expanding regions, when they have still a compa- 
rable size, but will be driven towards the value of the 
fastest expanding domains in the late-time limit. This 
growth in the average expansion rate corresponds to an 
acceleration of the volume scale factor. 

A convenient way to illustrate the possible departure 
of the solutions of the equations for the average scale 
factor, from the Friedmann solution, is the phase space 
analysis of scaling solutions provided in [7J : the standard 
Einstein-de Sitter model is unstable, a saddle point in 
the enlarged class of solutions for averaged inhomoge- 
neous models. Perturbations of the homogeneous state 
in the matter-dominated era will therefore drive the av- 
erage properties of the Universe away from the standard 
model in the direction of an accelerated, expanding and 
almost-isotropic state. This is also the region in which 
the Qx) term is positive. Therefore the instability of the 
Friedmann model leads naturally to accelerated expan- 
sion, if the phase space is traversed rapidly enough. A 
more elaborated phase space analysis may be found in 
L21J. 

The important changes to the Friedmann model that 
emerge when passing to explicit averages, may be sum- 
marized in the following list of generalizing concepts: 

1. The homogeneous and isotropic solution of Ein- 
stein's equations is replaced by explicit averages of 
the equations of General Relativity. 

2. The description is background-free, and inhomo- 
geneities do no longer average out on a background. 

3. The averaged model can be regarded as a general- 
ized background that is generically interacting with 
structures in the spatial hyper surfaces. 

4. The full Riemannian curvature degree of freedom 
is restored and the geometry no longer is restricted 
to a constant-curvature space. 

A more detailed review of our current understanding of 
the description of the average evolution may be found in 

Partitioning models 

To build a specific model using the above framework 
there have been several attempts. As the equations for 
ax) {t) are not closed, one has to impose, like in the 
Friedmann case, an equation of state for the fluid. The 
problem is here, that the fluid composed of backreac- 
tion Qx and average curvature (71) is only an effective 
one. Therefore it is not clear which equation of state 
one should choose, since it is dynamically determined by 
the inhomogeneities that are averaged over. In [20j for 
example, the equation of state of a Chaplygin gas has 
been used. Another approach in [7J has been to take a 



constant equation of state that is popular in quintessence 
models. This leads to simple scaling solutions for the ax 
dependence of Qx and {Tl)j^. Those will be generalized 
here by the partitioning model, described in detail in [23]. 



Model construction 

As the name suggests, our model implements a general 
partitioning of a domain of homogeneity V into subdo- 
mains [8j. To have a physical intuition about the evolu- 
tion of the subdomains, a reasonable choice is to partition 
V into over-dense M- and under-dense £^-regions. M 
and 8 regions will also obey the average equations (2) 
and there are consistency conditions, resulting from the 
split of the V equations into M and £ equations, that 
link the evolution of P, M and £. The reason for the 
split is that it offers the possibility to replace the less ac- 
cessible backreaction parameter Qx by a quantity that 
illustrates more directly the departure from homogeneity, 
namely the volume fraction of the over-dense regions \m • 
As explained above, \m is expected to decrease during 
the evolution and it is this decrease that drives the accel- 
eration. The physical motivation to split into over- and 
under-dense regions is that from the structure of the cos- 
mic web one may expect 8 regions, which are mainly com- 
posed of voids, to be more spherical than M. regions. On 
f , the expansion fluctuations should therefore be larger 
than the shear fluctuations and so Qg should be positive. 
The shear fluctuation-dominated M. regions should have 
a negative Qm- This increases the difference between 
the faster expanding £ and the decelerating M regions 
even further and therefore magnifies the expansion fluc- 
tuations on V that drive acceleration via a positive Qx- 
The fact that Qm and Qs are nonzero is the main dif- 
ference to a similar model of Wiltshire [241 [25], 

The generalization of the single scaling laws mentioned 
above is achieved by imposing the scaling on AA. and on 
£. In [14J the authors showed that Qx and {lZ)j^ may 
be expressed in a Laurent series starting at a'^ and a'^ 
respectively. This behavior breaks down if the fluctua- 
tions with respect to the mean density become of order 
one. This happens later on M. as well as on because 
on these regions the mean values lie above, respectively 
below the global mean, and therefore the fluctuations on 
Ai with respect to this over-dense mean are smaller than 
the variation between the peaks on Ad and the troughs 
on E. Therefore, the -scaling for Qx on M and 8 is 
expected to hold true even if the V regions already de- 
part from this perturbatively determined behavior. The 
partitioning model is therefore the first step of a general- 
ization to a physically justified nonlinear behavior of Qx 
and (JZ)j^ on the global domain V. 

Using this ansatz for the Qx)-evolution on M. and ^, 
and exploiting the consistency conditions for the parti- 
tioning, it is possible to arrive at a model that depends 
only on three parameters: the Hubble rate today Hxq^ 
the matter density today and the volume fraction of 
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FIG. 2: Evolution of the backreaction parameter Qt> in the 
course of cosmic time. From a value of of around 10~^ 
at a redshift of z ^ 1000, Qt) first decreases like and 
becomes approximately constant for a long period of the evo- 
lution when Xm gets important. 



the over-dense regions today Xmo- this parametriza- 
tion i^x^o o^ly sets the time scale, so that we may fix the 
evolution with and Xmo • Assuming for the con- 
cordance value of 0.27, the model shows that the more 
structure there is, indicated by a low value of Xmo^ 
higher the acceleration on the domain V will be. 

To analyze what the order of magnitude of Xm is to- 
day, an N-body simulation was studied. Smoothing the 
point distribution on 5h~^Mpc and applying a simple 
number count for the determination of Xmq^ a value of 
0.09 was obtained. An analysis by a Voronoi tessellation 
gave about 0.02. To obtain a definite value, this analy- 
sis clearly has to be improved. One would have to use a 
proper SPH smoothing and trace the Lagrangian regions, 
that are fixed in the initial conditions, until today. But in 
any case, Xmq seems to be in the region below 0.1. Inter- 
estingly enough, this is the region for which leads to 
a nearly constant Qx) on P, shown in Fig. 2. Hence, for 
low values of Xmq, Qv acts like a cosmological constant. 



Observational strategies 

This may also be seen by a fit of the model to luminos- 
ity distances of supernovae (SN). To convert the evolution 
of the average scale factor ax) into luminosity distances, 
we use the result of Rasanen [18j, who investigated the 
propagation of light in inhomogeneous universe models 
and provided a formula linking the observed luminosity 
distance to the volume scale factor ax)- The resulting 
probability contours in the parameter space Q'^^-Xmo in 
Fig. 3 show, that indeed the region around an of 0.3 
and a below 0.1 is favored by the data. 

To further test the viability of the model, it will soon 
be compared with more observational data. But as it is 
probably possible to fit also these, as indicated by the 
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FIG. 3: Probability contours of a fit [Martin Kunz priv. 
comm.] of the partitioning model to the Union2 SN data 
set. The model gives a comparably reasonable fit as a simple 
ACDM model. 

success of [13j, we need a quantity that will definitively 
allow to find out whether this model or one with a real 
cosmological constant fits better. To decide that, one 
may use a quantity introduced by Clarkson [9j, namely 
the C-function 

C{z) = l^H^ {DD" - D'^) + HH'DD' . (3) 

It consists of derivatives of the Hubble rate H {z) and 
the angular diameter distance D{z)^ and is constructed 
such that for every FRW model it is exactly for any 
z. For the partitioning model this is generically not the 
case and for several choices of parameters the difference 
is shown in Fig. 4 using the template metric of [13]. Un- 
fortunately, the function C (z) is too complicated to eval- 
uate it using present-day data, but as shown in [13], the 
EUCLID mission may be able to derive its values. 

Conclusion 

Routing the accelerated expansion back to inhomo- 
geneities furnishes an interesting possibility to avoid 
problems with a cosmological constant, such as the co- 
incidence problem and would give the sources of acceler- 
ation a physical meaning. If perturbations are analyzed 
on a fixed FRW background it is clear that the effect is 
way to small to give rise to an acceleration on the scale 
of the Hubble volume ([1, 2j). There are, however, semi- 
realistic non-perturbative models [17J that show that a 
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FIG. 4: Plot of Clarkson's C-function for several models dis- 
cussed in [23 . The model presented here is shown as the 
dotted line. The dashed line is its non-pert urbative general- 
ization explained in [23 , and the solid line is a single scaling 
model. For every Friedmann model C (z) = 0, which allows 
to distinguish averaged models from FRW models. 



considerable effect is not excluded. Furthermore it has 
been shown why perturbative models are not able to give 
a definite answer on the magnitude of the effect [ll], [5]. 
It seems that the question of how the growth of struc- 
ture influences the overall expansion of the Universe is 
still an open issue. The presented partitioning model can 
demonstrate that a ACDM-like expansion is possible in 
the context of these models, without the prior assump- 
tion of Qx) = const. ^ which emerges here more naturally. 
Furthermore it has been shown how this is related to the 
formation of structure described by the parameter X^. 
Finally, using Clarkson's C-function, we will have, at the 
latest with the data from the EUCLID mission, a handle 
on how to distinguish acceleration due to inhomogeneities 
from acceleration due to the presence of a strange fluid. 
All in all, these are promising prospects for the future. 
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